Formation of Liesegang patterns in the presence of an electric field 
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The effects of an external electric field on the formation of Liesegang patterns are investigated. 
The patterns are assumed to emerge from a phase separation process in the wake of a diffusive 
reaction front. The dynamics is described by a Cahn-Hilliard equation with a moving source term 
representing the reaction zone, and the electric field enters through its effects on the properties of 
the reaction zone. We employ our previous results [I. Bena, F. Coppex, M. Droz, and Z. Racz, J. 
Chem. Phys. 122, 024512 (2005)] on how the electric field changes both the motion of the front, 
as well as the amount of reaction product left behind the front, and our main conclusion is that the 
number of precipitation bands becomes finite in a finite electric field. The reason for the finiteness 
in case when the electric field drives the reagents towards the reaction zone is that the width of 
consecutive bands increases so that, beyond a distance £+, the precipitation is continuous (plug is 
formed). In case of an electric field of opposite polarity, the bands emerge in a finite interval £-, 
since the reaction product decreases with time and the conditions for phase separation cease to 
exist. We give estimates of £± in terms of measurable quantities and thus present an experimentally 
verifiable prediction of the " Cahn-Hilliard equation with a moving source" description of Liesegang 
phenomena. 

PACS numbers: 05.60.Gg, 64.60.Ht, 75.10.Jm, 72.25.-b 
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I. INTRODUCTION 

Precipitation patterns named after Liesegang Q, Q 
have been investigated for more than 100 years. Their 
discovery came from work on photoemulsions p|, and 
the interest in these patterns was sustained during all 
these years by recognizing that the underlying dynamics 
had connections to rather diverse physical- (e.g. near- 
equilibrium crystal growth 0), geological- (formation 
of agates 4J), chemical- (pattern formation in reaction- 
diffusion systems §) and more exotic (e.g. aggregation 
of asphaltene in crude oil @) phenomena. From theoret- 
ical point of view, the Liesegang patterns enjoyed contin- 
uing attention since attempts at their explanations were 
testing the theories of precipitation processes 0, H, 
(for a recent overview see Ref.|l(j) and, furthermore, 
the phenomenon was considered as a highly-nontrivial 
exam ple of p attern formation in the wake of a moving 

front mm. 

The origin of recent attention to Liesegang patterns is 
the hope that the phenomena may be relevant in engi- 
neering of meso- and microscopic patterns 0, 0, E| . 
The novelty of the idea is that, in contrast to the "top- 
down" processing (removing material in order to create a 
structure), the Liesegang dynamics provides a "bottom- 
up" mechanism where the structure emerges from a bulk 
precipitation process 0]. Of course, many obstacles will 
have to be overcome before such a strategy succeeds, the 
main one being the problem of controlling the pattern 
generated by a reaction-diffusion process. It is known ex- 
perimentally, with the results formulated in the Matalon- 
Packter law ^(|, that some degree of control may be 
exercised through the appropriate choice of the concen- 
trations of the inner and outer electrolytes participating 



in the process. It is also known, but much less under- 
stood, that the gel strongly influences the resulting pat- 
terns | l7| . Furthermore, recent experiments |15j indicate 
that the shape of the gel may be used in designing ap- 
propriate geometry in the precipitation patterns. 

The methods of control described above are somewhat 
rigid since the parameters cannot be changed during the 
process while, ideally, one requires an easily tuned, flex- 
ible external field for control. In principle, the electric 
field provides such an external control and, indeed, there 
is experimental evidence [lllliiSplllllillilillll 
HEHEHMHOl that an electric field significantly alters the 
emerging pattern (see Fig^ for a schematic experimen- 
tal setup). Unfortunately, these experiments give only a 
qualitative picture about the effects of the electric field, 
and the theory is even less developed pol I27L I2ll2ll3ll| . 

Our aim with this paper is to improve the theoreti- 
cal description of electric field effects and to bring it up 
to the level of quantitative predictions. We have shown 
recently l32| th at there are problems with previous at- 
tempts |20l l27l l28l 12^. l3l| which incorporate the electric 
field by assuming that it results in the drift of the reacting 
ions as well as in the drift of the reaction zone. Treating 
the background ions properly by using the electroneutral- 
ity condition, we found that the main effect of the field 
is that the amount of reaction product left behind the 
reaction zone increases (decreases) linearly in space de- 
pending whether the field drives the reacting ions towards 
(away) from the reaction zone. Building on these results, 
we shall show below using the "Cahn-Hilliard equation 
with moving source" model [33l l34l] , how the field affects 
the precipitation pattern itself. 

The choice of the model must be explained since a 
century of research did not lead to a generally accepted 
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theory of this pattern forming process. The reason is 
perhaps the complexity of the interplay between the mo- 
tion of the reaction front and the precipitation dynamics 
of the reaction product (and the intermediate reaction 
steps that may also be present), thus preventing the cre- 
ation of a single model encompassing all the possibilities. 
The approach we consider 33, 34] simplifies the situation 
by first treating the kinetics of reaction and the motion 
of the reaction front |3^, |3(| |33|. Then, the reaction 
product generated by the front is inserted as a source 
in the phase separati on p rocess described by the Cahn- 
Hilliard equation 133, 139, 0, E] ■ This theory has been 
shown [33l |3 0) EI El to generate Liesegan g pa tterns 
which satisfy the time- and spacing-laws [4a, |4y] (pat- 
terns emerging from most other theories do the same), 
the Matalon-Packter law 0] (significantly fewer theories 
can produce such patterns), and the width law p^lrjUIsf 
(no other theory yields this law since none of them has 
an underlying thermodynamics for setting the values of 
steady-state concentrations). This theory is also distinct 
from previous approaches in that it has a minimum num- 
ber of parameters which can be related to experimentally 
measurable quantities and thus the theory can give quan- 
titative predictions 42]. We should also note that, de- 
pending on the motion of the front, the phase separation 
process may occur at the position of the front or well be- 
hind the front [34| thus the model can also describe the 
limit of small imposed gradients |43 • 

In order to summarize the results obtained from the 
above theory for the case of external field present, let us 
first describe the setup (FigQ more precisely. A chem- 
ical reagent B = (B + ,B~), called inner electrolyte, is 
dissolved in a gel matrix inside a part of the small cen- 
tral cylinder (from point "0" to the right on the m-axis). 
A second reactant A = (A + ,A~) (outer electrolyte), of 
much higher concentration, is brought in contact with 
the gel (at point "0" on the x-axis). The two reservoirs 
of electrolytes A and B are providing constant concen- 
trations do and &o (o-o ^> bo) of the reagents at the two 
ends of the central cylinder. The outer electrolyte A dif- 
fuses into the gel and reacts (A~ + B + — > C) with the 
inner electrolyte B. The reaction front moves to the right 
and, under appropriate conditions, the reaction product 
C precipitates, and one observes the emergence of bands 
of precipitate perpendicular to the direction of motion 
of the front. The reservoirs of electrolytes A and B are 
kept at a constant potential difference U = Vb — Va, 
which corresponds to an average applied field intensity 
E = — U/L inside the central cylinder of length L. Let 
us note that throughout this paper we shall call forward 
electric field a field that drives the reacting ions towards 
the reaction zone; in our setup this corresponds to a pos- 
itive tension U > 0, or to E < 0. On the contrary, a field 
that works against the reacting ions reaching the reaction 
zone will be referred to as a reverse electric field, which 
in our setup corresponds to a negative tension U < 0, or 
to E > 0. 

Without the external field (U = 0), the positions of 
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FIG. 1: Schematic representation of the system under study. 
The inner electrolyte B(B + , B~) is dissolved in a gel matrix 
inside a part of the small central cylinder (from point "0" to 
the left on the s-axis). The outer electrolyte ^4(^4+, A~), of 
much higher concentration, is brought in contact with the gel 
(at point "0" on the x-axis). The two reservoirs of electrolytes 
A(A + , A~) and B(B + , B~) assure constant concentrations ao 
and 60 («o ^ bo) of the ions at the two ends of the central 
cylinder. The outer electrolyte ions diffuse into the gel, where 
the reaction A" + B + — > C takes place. (The background 
ions A + and B~ do not react.) The reaction front moves 
without convection towards the right, and precipitation bands 

- i.e., the alternation of high-density- C regions (shaded areas) 
and low-density-C regions - emerge inside the cylinder in the 
wake of this moving reaction front. A tension U = Vb — Va is 
applied between the ends of the central cylinder of length L, 
which corresponds to an average applied field intensity E — 

— U/L. A positive tension U > (corresponding to E < 0) 
drives the reacting ions towards the reaction zone (forward 
electrc field). A negative tension U < (corresponding to 
E > 0) works against the reacting ions reaching the reaction 
zone (reverse electric field). 



the bands (x n , measured from the initial contact of the 
reagents) form a geometric series, x n ~ (1 +p) n , where 
1 + p or p > is called the spacing coefficient. For small 
fields (\E\ = \U\/L < 2 V/m, for which one has a 'suf- 
ficient' number of bands, i.e., > 20), the band spacing 
can still be described as geometric series and our first 
result pertains to the dependence of the effective spacing 
coefficient on the applied field, p = p(U/L). We call it ef- 
fective spacing coefficient because these geometric series 
are finite as evidenced by the results for higher fields. 
The reason for the finite number of bands for electric 
field that drives the reagents towards the reaction zone 
(called forward field in the following) is that the width 
of consecutive bands increases faster than the distance 
between the bands. Thus one finds that the precipita- 
tion is continuous (plug is formed) beyond a distance 1+ . 
The number of band is also finite and they appear in an 
interval [0, tJ\ for the case of a field of reverse polarity. 
The reason for finiteness, however, is different. In case of 
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a field working against the reacting ions reaching the re- 
action zone, the amount of reaction product generated in 
the reaction zone decreases with time and the conditions 
for phase separation cease to exist. The quantities and 
£- are easily accessible in experiments, and our main re- 
sult (apart from qualitative observations) is that we give 
an estimate of them in terms of measurable quantities. 
Thus we provide a way of designing the spatial range of 
the emerging pattern. In view of the competing theories 
of Liesegang phenomena, this also gives yet another way 
of discovering which is the correct description. 

Below we present the details in the following order. 
First, the Cahn-Hilliard equation and the changes in the 
source term in the presence of an electric field are dis- 
cussed (Sec. II). Next (Sec. Ill), the results of numeri- 
cal simulations and qualitative arguments concerning the 
characteristics of the Liesegang patterns for a forward 
applied field are presented. The case of the reverse po- 
larity is described in Sec. IV. The comparison with ex- 
periments is discussed in Sec. V, while the conclusions 
and perspectives are given in Sec. VI. Finally, a discus- 
sion on the choice of the parameters of our theoretical 
model, as inferred from experimental data, can be found 
in the Appendix. 



II. PHASE SEPARATION AND DYNAMICS OF 
THE REACTION PRODUCT C IN THE 
PRESENCE OF AN EXTERNAL ELECTRIC 
FIELD 

There is experimental evidence that the Liesegang pat- 
tern is an alternation of high-density Ch and low-density 
ci regions of the precipitate C. Thus, the reaction prod- 
uct phase-separates behind the reaction front, and since 
this takes place on a macroscopic scale, its dynamics 
can be represented through an extension of the Cahn- 
Hilliard equation (in other contextSjthis is the equation 
of model B of critical dynamics |41|). This equation, 
however, requires the knowledge of the free energy den- 
sity T of the system; note that for a homogeneous system, 
T has to present two minima in order to accommodate 
the two equilibrium states of high Ch and low ci densi- 
ties. The simplest form of T having this property and 
containing the minimum number of parameters (which, 
moreover, guarantees the stability of the system against 
short-wavelength fluctuations) is the Ginzburg-Landau 
free energy density 



J-[m] 



-em 2 + -jm 4 + -<j(\7m) 2 . (1) 
2 4 2 



Here we introduced the 'reduced' concentration 
c{x,t) - (c h + ci)/2 



m(x, t) 



(c fc - cj)/2 



(2) 



that varies between —1 and +1 when the concentration 
of the product C varies between ci and c/,. Note the un- 
derlying assumption of the one-dimensional character of 



the system, i.e., the fact that all the relevant parameters 
are only ^-dependent, that is well-justified for the exper- 
imental setup in Fig. ^ The parameters e, 7, and a are 
system and temperature dependent, with (i) e > en- 
suring that the system is in the phase separating regime 
(i.e., the temperature is smaller than the critical one); 
(ii) 7 = e in order to ensure that the minima of the free 
energy correspond to m = ±1 (the homogeneous high 
and low density phases of C) ; (iii) a > in order to pro- 
vide stability against short-wavelength inhomogeneities. 
Figure |21 offers a schematic representation of the homo- 
geneous part of the free energy J 7 as a function of to, 
with the different stability regions. Because the C-s are 
assumed to be neutral particles, the free energy density 
(i.e., the parameters e = 7, a, and A) and the corre- 
sponding stability diagram are not modified by an exter- 
nal electric field. 




FIG. 2: The homogeneous part of the free energy density T as 
a function of the reduced density m. The phase separation is 
an activated process in the metastable regimes m s — 1/ \/3 < 
\m\ < m e = 1, while it goes by spinodal decomposition in the 
linear instability region \m\ < m s = 



The dynamics of C is thus driven by the free-energy 
density l|T]). but, in addition, there is a continuous cre- 
ation of C by the moving reaction front, with a certain 
space and time dependent source density. One is thus led 
to write down |33| the following phenomenological evolu- 
tion equation for the reduced concentration field m(x, t): 

= -AA(em-7TO 3 + crATO,) + S' m (a;,i), (3) 
where A is a kinetic coefficient and S m (x,t) is the source 
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density. Equation 10 has to be solved with the homoge- 
neous initial condition m(x, t = 0) = — (c^ + c{)/(ch — c/) 
that corresponds to the absence of C inside the system 
before the beginning of the reaction. Note that the above 
"Cahn-Hilliard equation with a source" should also con- 
tain two noise terms. One of them is the thermal noise, 
while the other one originates in the chemical reaction 
that creates the source term. However, as discussed 
in [3(| , the noise in A+B — > C-type reaction fronts can be 
neglected in dimensions d > 2, while neglecting the ther- 
mal noise term means that an effective zero-temperature 
process is considered, and that the phase separation takes 
place only through a spinodal decomposition mechanism. 
This approximation is supported by the experimentally- 
known fact of the very long life of the formed patterns, 
which amounts to a very low 'effective temperature' of 
the system. The theory could be refined by including the 
thermal noise, since then the nucleation and growth pro- 
cesses would be also captured. The role of noise has been 
investigated for the fieldless case in [34|], and its effects 
in the presence of an electric field will be the subject of 
a forthcoming paper. Here we shall remain within the 
deterministic framework corresponding to Eq. 

Let us now concentrate on the source term S m {x, t) 
in Eq. J3J. As already mentioned, this term models the 
production of C by the moving reaction front, and it is 
influenced by the presence of an external electric field. 
The effect of the electric field has been studied in de- 
tail in [82l |. and we summarize below the main results 
for the range of parameters that are relevant for typical 
experimental situations. One has to realize first that the 
reagents A and B are electrolytes which dissociate, 



A -» A+ + A~ 



B — > B 



B~ 



and the basic reaction process is 



A' 



B^ 



C. 



(4) 



(5) 



while the 'background' ions A + and B~ are not react- 
ing. The modelling of the system in |32| was based on 
several simplifying assumptions: (i) the one-dimensional 
character of the system (i.e., all the relevant parameters 
are only ^-dependent); (ii) the complete dissociation of 
the electrolytes A and B into their respective ions (that 
allows to eliminate the dynamics of the neutral As and 
B-s from our description); (iii) infinite reaction rate and 
irreversibility of the basic reaction A~ +B + — > C. This is 
justified by the fact that the characteristic reaction time- 
scale is much smaller than any time-scale connected with 
diffusion and precipitation (pattern formation) , and leads 
to a point-like reaction zone; (iv) the electroneutrality 
approximation (the local charge density is zero on space 
scales that are relevant to pattern formation), whose ap- 
plicability for the systems under study was discussed in 
detail in |37| : (v) we considered monovalent ions; (vi) 
finally, we assumed equal diffusion coefficients D of the 
ions. As mentioned in [32| . any of these assumptions may 
be relaxed without generating major changes of the con- 
clusions of our study. The reaction-diffusion equations 



for the concentration profiles of the ions can be solved 
numerically (with boundary conditions that correspond 
to the presence of the two reservoirs of ions - of con- 
centrations ao and bo, respectively - at the ends of the 
reaction cylinder). 

A first result refers to the motion of the point-like re- 
action zone; it is, to a good approximation, a diffusive 
motion, 



x f (t) = V2Dji, 



(G) 



with a diffusion coefficient Df that is practically unaf- 
fected by the field intensity, i.e., it is given by its fieldless 
expression 




(ao/M 



1 



Df_ 

2D j {a Q /bo) + 1 ' 



(7) 



where D is the common diffusion coefficient of the ions. 
Note that for a reverse field there is also a small drift 
component in the motion of the front. However, for the 
range of reverse fields and observation times we are con- 
sidering, this drift component can be neglected. 

The effect of the field is more significant in the pro- 
duction of C. As known [ljj, in the absence of the field 
the concentration of C-s left behind the front is a con- 
stant and its value cq is determined by the initial con- 
centrations of the ions ao and bo, and by their diffusion 
coefficients. In the particular case of equal diffusion co- 
efficients D of the ions, its value is given by: 



co^a a KJ2D/Df 



(8) 



where K = (1 + 6 /ao)(2 v / 7r)" 1 exp(-D f /2D), and the 
diffusion coefficient Df of the front is given by Eq. J7J. 
For an infinite reaction rate, the corresponding source 
density for the production of C (i.e., the density of pro- 
duction of C per unit time) is a 8 peak localized at the 
instantaneous position of the front Xf(t), 



Sjnix, t) 



a KVD 5(x-x f (t)) 



{c h ~c{)/2 



Vt 



(9) 



whose amplitude decays in time as ~ ^-/yi. 

Consider now an external electric field applied to the 
system. For both polarities of the field, and for relatively 
small values of the field intensity (the so-called linear 
regime at which we shall limit our study in this paper), 
there is a linear variation of the concentration of C-s with 
x, with a slope that is proportional to the applied tension, 



c(x)=cq[1 + (7]U/L)x], 



(10) 



as illustrated by Fig. |21 The parameter -q depends on 
the values of the other parameters (i.e., a , b , and D), 
and for the set of parameters in Fig. |3| one infers rj w 5 
V -1 . We should mention that the validity of this linear 
regime is wider for positive tensions (e.g., it may go up to 
U/L w 10 V/m for the system considered in Fig. and 
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less extended for negative tensions (e.g., up to U/L « —2 
V/m for the system in Fig.[3J); beyond these limits there 
is a relative error larger than 10 % in approximating c(x) 
by a linear profile with the initial slope. 




FIG. 3: The density of the reaction product C left behind 
the reaction front for different values of the electric field U/ L 
applied to the system. The observation time is t = 10 days. 
The values of the other parameters are: ao = 10 M, bo =0.1 
M, and D = 10~ 9 m 2 /s, that lead to c = 0.1145 M. 



The above result on the spatial dependence IjlUI) of the 
concentration of the reaction product is incorporated into 
the source term through the following modification of its 
amplitude: 



a KVD [l + (r)U/L)x]8{x-Xf(t)) 



(<* - c t )/2 



(11) 

As discussed in the Introduction, in the absence of an 
electric field (U = 0) this spinodal decomposition sce- 
nario reproduces, in a simple and coherent way, all the 
generic laws of Liesegang patterns. Moreover, it contains 
very few parameters, which can be inferred from exper- 
imental data H2, and thus has a predictive power. We 
expect to recover these qualities in the presence of an 
applied external electric field, as well. 

The solution to the Cahn-Hilliard equation © with the 
source (|llfl . i.e., the profile of the reduced concentration 
m(x, t) is obtained numerically. As already mentioned 
above and discussed in detail in |32, we decided to fo- 
cus our analysis on the situations that are experimentally 
relevant and, in particular, the choice of parameters in- 
tends to mimic real experimental situations. Namely, we 
considered concentrations of the reagents ao and bo in 
the range 10 -2 — 10 M, length L of the system of some 



tenths centimetres, and tensions U applied between sys- 
tem's edges such that we are in the linear regime of the 
production of C, i.e., U/L varies between —2 and +10 
V/m. The common diffusion coefficient of the ions was 
chosen as D — 10 -9 m 2 /s for all the calculations. The pa- 
rameters A, e = 7, and a of the free-energy-driven part of 
the dynamics of C can be inferred from the experimental 
data as explained in the Appendix. Finally, the pattern 
formation was followed for a period of the order of ten 
days of experimental observation time. In the following 
section we present the results that were obtained for a 
polarity U > of the applied field (the forward field), 
that favours the reaction, i.e., that drives the A~ and 
B + ions towards the reaction zone. Section IV will be 
devoted to the U < (the reverse field) case study. 



III. PATTERN CHARACTERISTICS FOR A 
FORWARD APPLIED FIELD 

A first generic feature of the patterns formed in the 
presence of the forward field (U > in our setup) is 
that the number of bands is finite, i.e., band formation 
stops at a certain moment through the appearance of a 
continuous precipitation region (a 'plug' of high-density 
precipitate). The higher the tension, the earlier this plug 
forms, see Fig. 0] for an illustration. 

One can make a rough estimate of the distance 1+ of 
the onset of the plug through the following reasoning: 
the plug forms when the width of the n-th high-density 
band, w n , becomes equal to the distance (x n +i — x n ) 
between the n-th and the (n + 1) high-density bands. 
If we estimate that the amount of C produced by the 
front between x n and x n+ i goes entirely in the n-th high- 
density band, i.e., 

c [l+(x n +x n+ i)/2](x n+1 -x„) = c h w n = c h (x n+1 -x„) , 

(12) 

then we obtain for the distance £ + = (x n +x n+ i)/2 « x n : 



L Ch 
rjU 



c 



CO 



(13) 



Figure [S] offers a comparison between the results of the 
numerical simulations and this theoretical estimation of 
£ + for different values of the forward field U/L > 0. Note 
that the experimental measurement of £ + allows us to 
estimate the value of Ch (provided that Co is known). 

A second generic feature of the pattern is that the dis- 
tance between two successive bands diminishes as com- 
pared to the fieldless case, and this effect is increasing 
with increasing forward applied field U/L. This can be 
easily understood through a simple qualitative argument. 
In the presence of the forward field, the reaction front 
leaves behind a larger quantity of C than in the absence 
of the field. Thus, after the formation of a band, the 
re-establishment of the phase-separation instability con- 
ditions takes place sooner in the presence of the field, 
resulting in a higher spatial density of bands in the sys- 
tem. With a good approximation, the positions of the 
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FIG. 4: The profile of the reduced concentration m(x) for 
different values of the forward applied field U/L > 0. The 
snapshots are taken at t = 14 days. The dashed lines repre- 
sent the position of the reaction front at this time. The values 
of the other parameters are Co = 0.1145 M, ci = 0.0045 M, 
c h = 0.3645 M, D = 10~ 9 m 2 /s, D f = 5.43 ■ 10" 9 m 2 /s, and 
rj — 5 V - . One notices the decrease in the band spacing 
with increasing tension, as well as the appearance of the plug 
- earlier for larger tensions. 



FIG. 5: The value of the distance £+ of the plug-setting-in as a 
function of the forward applied field U/L > 0. The diamonds 
represent the numerically estimated values, while the dashed 
fine corresponds to the rough theoretical estimate, Ea. ll3> . 
The values of the parameters of the system are the same as 
for Fig. 01 The size of the symbols is representative for the 
estimated error bars. 



bands x n still form a geometric series as in the fieldless 
case, and one can define an experimentally measurable 
'effective' spacing-law parameter p: 

x n ~{l+p) n . (14) 

Thus, as illustrated in Fig. [SI p is a decreasing function 
of the forward applied field U/L. Due to the decrease in 
the number of bands with increasing field intensity our 
measures of p were restricted to a rather narrow interval 
of field intensities around U/L = 0. 

Another remarkable feature of the pattern is its 'insta- 
bility'. Indeed, as illustrated by Fig. [7] due to the high 
densities of C in the plug region (which thus becomes 
unstable), there is a 'flow-back' of C from the plug to- 
wards the regions where the patterns had formed previ- 
ously. Therefore, this may cause a gradual disappearance 
of some of the already formed bands. This effect is not 
present in the fieldless case. 

One has to realize that the structure of the patterns 
depends essentially on the position of the concentration 
of C-s left behind by the reaction front with respect to the 
stability domains associated with the free energy density 
T . This remark holds both in the absence and in the 
presence of an applied field. If, for example, cq is below 
the spinodal value, no pattern will form in the absence of 



the field. However, in the presence of a forward field, the 
concentration of C-s left behind the front is continuously 
increasing, and thus, at a certain moment, it will cross 
the spinodal curve, and therefore the pattern formation 
mechanism will be 'turned-on', i.e., patterns will start to 
form after a certain time and then continue up to the 
eventual plug formation discussed above, as illustrated 
in Fig. Hi. 

We can also comment here on the width law 0> 03 
Egj . As shown in (4^, in the fieldless case the derivation 
of this law, w n ~ x", with a = 1, is straightforward. Let 
us apply the same reasoning as in in the presence of 
an electric field. One combines the facts that (i) the re- 
action front leaves behind a density c(x) of C-s, Eq. 1(10)1: 
(ii) the C-s segregate into low (ch) and high (ch) density 
bands; (iii) the number of C-s is conserved in this segre- 
gation process. The equation expressing the conservation 
of C-s can be written as 

c [1 + (rjjU/L){x n + x„ + i)/2] (x„ +x -x n ) = c h w n +ci(x n+ i 
which together with x„+i = (1 + p)x n leads to 

wn = Mcp^Q) ^ + c p(l+p/2) V U/L xl 
Ch — ci c h - Cl 

It is however difficult to make a clear-cut statement on 
this point. Indeed, for field intensities for which the sec- 
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FIG. 6: The value of the effective spacing coefficient p as a 
function of the applied field U/L. The values of the parame- 
ters of the system are the same as for Fig. 2] The size of the 
symbols is representative for the estimated error bars. 



FIG. 7: The profile of the reduced concentration m(x) at dif- 
ferent times for a fixed forward applied field U/L — 6 V/m. 
The dashed lines represent the position of the reaction front at 
the corresponding times. One notices the progressive 'flow- 
back' of the C-s from plug zone, i.e., the backward disap- 
pearance of the pattern. The values of the parameters of the 
system are the same as for Fig. 0] 



ond term in l|16l) may start to play a role, the total num- 
ber of bands is so small that no reliable conclusion can be 
drawn about the 'systematic' behaviour of their width. 
In case one attempts to fit w n ~ x", one will infer an 
effective a which is between 1 and 2 and which increases 
with the field. 



IV. PATTERN CHARACTERISTICS FOR A 
REVERSE APPLIED FIELD 

Let us consider now the case of the reverse polarity 
of the electric field (i.e., the case when it drives the re- 
acting ions away from the reaction zone, U < in our 
setup) . Again, a generic characteristic of the pattern is 
the finite number of bands. This can be easily under- 
stood in connection with the continuous decrease of the 
concentration of C-s left behind the reaction front, up 
to a point when the phase-separation conditions are no 
longer fulfilled. As illustrated in Fig. |SJ there are less 
and less bands for larger and larger field intensities. But 
contrary to the forward field case, these bands are 'sta- 
ble' i.e., nothing analogous to the 'flow-back' process in 
this case. Moreover, the last formed band collects pro- 
gressively all the C-s left behind the front, i.e., its width 
increases slowly with time. 

There exists thus a maximum distance I- of the spatial 
extension of the pattern, which can be estimated experi- 
mentally as the right edge of the right-most high-density 
band. A rough theoretical estimation of this length is 



given through the condition that the concentration of 
the C product at this point equals the lower limit of the 
spinodal decomposition domain, i.e., 

from which 

e _ = _±_ (co-c t )-(l-l/V3)( Ch - Ct )/2 
r)\U\ c 

A comparison of this result with the numerical simu- 
lations is given in Fig. ^| for different values of the field 
intensity \U\/L. 

We can see that the expressions for the two lengths £ + 
and I- (Eqs. (|T3|l and ((T7|l % ) provide a way to evaluate ct 
and ci through experimental measurements of the length 
of the patterns. 

As expected, in the reverse field case, there is an in- 
crease of the spacing between the bands (as compared to 
the fieldless case), and this is illustrated by an increasing 
of the effective spacing coefficient p with increasing field 
intensity \U\/L, see Fig. for a numerical estimation. 

Finally, it is worth mentioning again how the appear- 
ance of the pattern is sensitive to the position of cq with 
respect to the stability domains associated with the free 
energy density T . If, e.g., cq is above this instability do- 
main, there is no pattern formation in the fieldless case, 
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FIG. 8: The profile of the reduced concentration m[x) at 
t = 14 days for two sets of values of the parameters: (a) Co = 
0.1145 M, D f = 5.43 • 1(T 9 m 2 /s, r\ = 5 V -1 , c h = 0.72 M, 
and c; = 0.120 M (i.e., Co is below the spinodal decomposition 
domain), and an applied electric field U/L — 10 V/m. One 
notices the initial absence of pattern, then a limited region of 
existence of the pattern and, in the end, a plug region starting 
to form, (b) c = 0.1145 M, D f = 5.43 • 10~ 9 m 2 /s, V = 5 
V" 1 , c h = 0.1045 M, and a = 0.045 M (i.e., c is above the 
spinodal decomposition domain), and an applied electric field 
U/L = —2 V/m. One notices the initial plug region, then a 
limited region of existence of the pattern, and, in the end, the 
disappearance of the pattern. The dashed line represents the 
position of the front at the time of the snapshot. 
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FIG. 9: The profile of the reduced concentration m(x) for 
different values of the reverse applied field U/L < 0. The 
snapshots are taken at t — 14 days. The dashed lines repre- 
sent the position of the reaction front at this time (except for 
the last panel, where it indicates the stopping of the reaction, 
S m — 0). One notices the increase in band spacing with de- 
creasing applied tension, as well as the rapid disappearance 
of the pattern. The values of the parameters of the system 
are the same as for Fig. 2] 



but just a uniform plug of C behind the front. How- 
ever, as illustrated in Fig.|HJj, in the presence of a reverse 
field, due to the continuous decrease in the concentration 
of C, the instability domain may be reached from above, 
and a pattern starts to form after an initial plug region, 
and eventually, as discussed above, stops after some time, 
when the instability domain is left for too small concen- 
trations of C . 



V. COMPARISON WITH EXPERIMENTS 

Experiments on Liesegang patterns in the presence of 
an electric field have been going on quite a while since 
the pioneering works of Happel et al. u3 and Kisch , 
see EH 112, Hp, SLUM EH, E3, E^Hp, [33, p3|. As al- 
ready mentioned in |32j| . with the exception of plll28Ll2^| . 
these experiments were carried on for a single polarity of 
the electric field (i.e., either a 'forward', or a 'reverse' 
field according to our terminology). [5l| It is therefore 



not surprising that each of them 'covers' only a part of 
the situations that are predicted by our theory, which 
thus has the merit to include all these elements in an 
unifying, coherent, and simple frame. We shall present 
below a brief qualitative summary of these experimental 
results. 

Concerning the motion of the front, it is experimentally 
found to be diffusive in [2l|, (2^ (with a slight decrease 
in the diffusion coefficient with increasing field intensity) , 
and (for a 2D geometry) , in agreement with our the- 
oretical predictions |33 | for a forward field. Other exper- 
iments [25l l2rl [27J , present a motion of the front with 
a small drift component (that increases with increasing 
field intensity |2l| ) . Occasionally, this drift component is 
given a theoretical justification on the basis of a reaction- 
diffusion model for the system with a s uperimp osed con- 
stant electric field intensity, see |27L l28l l29i 13 1| - that, as 
discussed in j^], is a somewhat unrealistic assumption. 
Here again our model is able to capture the appearance 
of this small drift component for the case of a reverse 
electric field. 

Several experimentally-studied modifications of the 
characteristics of the Liesegang patterns in the presence 
of an electric field are quoted in the literature, and 
sometimes they seem contradictory: 
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FIG. 10: The value of the distance £- of disappearance of 
the pattern as a function of the absolute value of the reverse 
applied field \U\/L. The diamonds represent the numerically 
estimated values, while the dashed line corresponds to the 
rough theoretical estimate, Eg. 11711 . The values of the pa- 
rameters of the system are the same as for Fig. [3] The size of 
the symbols is representative for the estimated error bars. 



field, see Section IV. 

(vii) Finally, a special remark on the results of [2^, |2t| , 
which sometimes match our results, sometimes are just 
the opposite of those obtained from our theory (as far as 
the polarity of the electric field is concerned) . In particu- 
lar, in |28| it was found that (a) the motion of the reaction 
front is diffusive for a reverse field and has a small drift 
component for a forward field (contrary to our theory); 
(b) the average spacing coefficient p decreases with in- 
creasing field (just like in our theory). In [2^| an attempt 
to fit the width law (see, e.g., [4^) in the presence of an 
electric field, w n ~ x™, leads to an exponent a that is de- 
creasing monotonically with the field - in opposition with 
the result of our discussion following Eq. ljTo|l . However 
it was argued [5(j that in his case the properties of the 
intermediate compounds are responsible for this 'anoma- 
lous' behaviour, and such effects are outside the range of 
our theoretical model. 

Most of the above examples show agreement with 
our theory at a qualitative level. It should be noted 
however that our model offers quantitative estimates of 
the changes as compared to the fieldless case, and thus 
is well-suited for direct comparison with the results of 
appropriately-designed experiments. 



VI. CONCLUSIONS AND PERSPECTIVES 



(i) A first striking observation is the appearance of 
an unifor m p recipitation for sufficiently high field- 
intensities |2Ct I2H I22I l23| and/or after a sufficiently long 
time |2l|. 

(ii) Also, there is an 'acceleration' of the band forma- 
tion [2]]], i.e., the time of appearance of the first band 
decreases with increasing field intensity p^j ]. 

(iii) A decrease of the spacing between successive bands 
(as compared to the fieldless case) with increasing field 
intensity is registered in [2(], |2^| . 

(iv) On the other hand, an increase of the spacing 
between successive bands is described in |2^, |2(| |27| . 
and [3(J (for a 2D experimental setup). 

(v) As reported in [23. 126L l3(i| , at a given time there are 
less bands formed for higher field intensities. 

(vi) The paper [2^ indicates a reduction of the initial 
"diffuse portion" with increasing field intensity, and [3(j 
also gives an example of the reduction of the initial 
"fuzzy zone" . 

One is now in position to recognize that all these fea- 
tures are recovered in a simple way within our theory. 
More precisely, (i)-(iii) and the first result in (vi) are ob- 
tained in our model in the case of a forward electric field, 
see Section III; while (iv), (v), and the second result in 
(vi) can be obtained theoretically for a reverse applied 



In this paper we examined the influence of an ap- 
plied electric field on the formation and characteristics 
of Liesegang patterns that appear in the wake of an 
A~ + B + — > C reaction front. The appearance of the 
pattern is related to the phase separation of the pre- 
cipitate C into a high and a low density phase. At a 
macroscopic, phenomenological level, the dynamics of C 
particles can be modelled using a Cahn-Hilliard equation 
supplemented with a space and time-dependent source 
term that describes the production of C-s by the reaction 
front. It was found that for both polarities of the applied 
field the pattern has a finite spatial extension. Moreover, 
measuring this spatial extension for various values of the 
applied tensions allows to infer the values of the concen- 
trations of C in the high and low density phases. The 
distance between the bands in the pattern is influenced 
by the intensity of the applied field, namely it decreases 
with increasing U/L for a forward field U > 0, and in- 
creases with increasing \U\/L for a reverse field U < 0. 

One has to underline that the proposed model contains 
a minimal number of parameters, and all of them can be 
inferred from experimental data. Thus the model, be- 
sides reproducing well the generic experimental laws of 
Liesegang patterns, has a clear predictive power that can 
be used to control experimental situations. The data we 
are presenting are the result of numerical simulations, 
and are based on a set of simplifying assumptions. It is 
expected that most of these assumptions do not affect 
our main conclusions. However, the details of the reac- 
tion process leading to the final precipitate C may be 
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more complex (e.g., implying several steps, intermediate 
compounds with different electric charges and diffusivi- 
ties, etc.), and the role of the electric field may also have 
unexpected features. The simplest extension of our work 
would be to consider a reaction scheme with one bivalent 
and two monovalent ions, a case that occurs in several ex- 
periments. Work on this type of systems is in progress. 

APPENDIX A: SETTING THE PARAMETERS 

Typical experimental situations correspond to concen- 
trations ao and 60 of the reagents of the order of 10~ 2 — 10 
M, with a ratio ao/60 ~ 10 — 100, and thus it is suit- 
able to choose the unit of concentration as n = 1 M. 
The diffusion coefficients of the reagents are of the or- 
der ~ 10 -9 m 2 /s, so that the length I and time r scales 
should be chosen such that £ 2 /t is of the same order of 
magnitude, 




Moreover, the experimental patterns have a total length 
of about 20 cm, and the time to produce such a pattern 
is of some tenth days - and this offers the order of mag- 
nitude of the diffusion coefficient of the reaction front, 
which is typically of the same order or an order of magni- 
tude larger than the diffusion coefficients of the reagents 
(depending on the ratio ao/b of the concentration of the 
reagents). The typical widths of the precipitation bands 
are of a few mm at the beginning, and approach ~ 1 cm 
at the end, and so are the distances between two succes- 
sive bands. From the visual observation of the beginning 
of the band formation it takes some tenth minutes for the 
band to be clearly seen, and some hundred minutes for its 
complete formation. Since the Cahn-Hilliard model has 
intrinsic length scale \foJe and time scale cr/(Ae 2 ), these 



have to be comparable to the typical length and time 
scales of the appearance of a single band. We are thus 
led to the following orders of magnitude for the length 
and time scales: 

t = yfcTJl ~ 10" 4 m, t = cr/(Ae 2 ) - 10s. 

There remains, however, the question of determining the 
concentrations of the two phases of the C precipitate, ct 
and c/. As explained in the main text, this can be done 
through measurements of the total spatial extent of the 
pattern in the presence of a forward, respectively reverse 
electric field. 

Here is the set of parameters that we used in most of 
the simulations presented above. For the concentrations 
of the reagents we take ao = 10 M, 60 = 0.1 M, and 
the common diffusion coefficient of the electrolyte ions is 
D = 10~ 9 m 2 /s. This leads, according to Eqs. J7J and 
(JHJ, to a diffusion coefficient of the front Df = 5.43 • 
10 -9 m 2 /s, and a concentration of the fieldless reaction 
product Co = 0.1145 M. For the above parameters the 
coefficient 77 in the expression (|10fl of c(x) is rj = 5 V -1 
. The parameters of the Cahn-Hilliard model are chosen 
such that I = y/cr/e = 10~ 4 m, and r = a j (Ae 2 ) = 40 s. 
Regarding the concentrations of the low and high density 
phases of C, we usually set them such that Co given above 
is in the instability domain of spinodal decomposition, 
c h = 0.3645 M, and cj = 0.0045 M. 
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